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ABSTRACT 



In order to determine the shape of the ionization curve in 
nuclear emulsions at values of T> 100 s we have measured the blob 
density of relativistic electron tracks. The pv of the pair pro- 
duced electrons used was determined by multiple scattering. A_ 
total of 84.45 centimeters of electron track was bl-<$b counted and 
scattered. From T = 100 to T - 5400 the data were combined into 
14 points each with a statistical uncertainty on blob density of 
less than 1%. These points indicate a level "plateau" and show 
no deviation from this plateau within our statistical accuracy. 
Pious were used to estimate the minimum of the ionization curve, 
and ratio of blob density plateau to blob density minimum is 
estimated to be 1.140 + .020. 
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1 . 



Introduction 



The classical theory of the energy loss of a charged particle 
passing through matter was developed by Bohr /!/ and later ex- 
tended quantum mechanically by Bethe 111 and Bloch /3 / . These 
theories predicted that the energy loss by ionization per unit 
path length was a function of velocity and charge. In particular 
the energy loss was to decrease approximately as X/v to a minimum 
at v/c - 0.96. After this minimum the energy loss should loga- 
rifchxmically increase indefinitely due to relativistic effects 
on the electric field of the charged particle. Swann/4 / predicted 5 
however, that the increase in energy loss would be limited, in 
condensed substances, by the polarization of the medium. Quanti- 
tative calculations by Fermi /5 / , treating the electrons in the 
medium as classic oscillators of one frequency, showed that the 
most probable energy-loss per unit path length should increase 
to and remain at some maximum value which depends on the density 
of the medium. This maximum value was therefore called the 
"Fermi Plateau 11 , and the polarization effect the "density effect." 
Fermi's theory has been extended and elaborated upon, using 
multi -frequency theory, by Wick /6/, Halpern and Hall /7/, 

Schdnberg /8 / , and Sternheimer /9 /; but these extensions show no 
significant difference at v— c. and also predict a "plateau." 
Although emulsion groups were eager to measure the relativistic 
and density effects in solid media, it was not until 1948 that 
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Berriman introduced photographic emulsions that were sensitive 
to singly charged particles of minimum ionization. 

The first attempts to measure the relativistic rise met with 

negative results. These early experiments are discussed by Price 

/I 0/ . However 5 in 1950 Pickup and Voyvodic /I 1 / reported the 

first experimental indications of a relativistic rise, Corson 

and Keck /12/, McDiarmid /13/ ? and Morrish /li / reported the 

existence of an ionization plateau j but none of these experiments 

carried their measurements down to the ionization minimum. Two 

more experiments published in 1952 showed a relativistic increase 

in grain density from minimum to plateau 5 but differed in the 

magnitude. Daniels et al / 1 5 / reported a rise of about 8%, 

while Shapiro and Stiller /16/ obtained an experimental ratio of 
+0.04 

1.12-0.03. Stiller and Shapiro Jill duplicated their experiment 
the following year using blob count instead of grain count as a 
measure of ionization and found virtually the same results-- 
a saturation of the curve at T>ioo, and a rise above minimum 
of about 1 Michaelis and Violet /18 / and Fleming and Lord 
/ 19 / also attempted to measure the ratio of ionization at plateau 
to ionization at minimum. In both experiments two points were 
plotted: one near minimum and one near plateau. These experi- 

ments have, underlying their design, the implicit assumption of 
a constant Fermi plateau. Experimenters were now trying to 
measure two parameters: at what value of 'Y' the ’’plateau 8 ' 



2 



begins, and how far above the ionization minimum the "plateau" 

lies. The results of Alexander and Johnston /20/ are quoted as 

a value of grain density plateau/grain density minimum of 

1.133 + .008, and Jongejans /21 / has stated in 1960: 

It seems that a remarkable agreement exists about 
the ratio plateau to minimum; we find 




1.129 4- .010 



This value compared with 1.133 + 0.008 of Alexander 

and 1.14 + 0.03 of Stiller. 

Recently an extensive study of the relativistic rise of grain 
density was carried out by Patrick and Barkas / 22/ s but thair main 
attention was focused on the rise from minimum and transition to 
plateau. 

Theoretical studies /23-26/ were continuing, but the major 
effort was devoted to finding the exact shape of the curve in the 
interval during the relativistic rise and transition to plateau. 

If the ionization curve could be plotted accurately in this range 
of OT , the identification of particles with velocities in this 
range could be made. The rate of energy loss per unit path could 
be measured; and together with a measurement of velocity, would 
give the particle identity even on tracks which passed from the 
emulsions. It is not surprising that the plateau region was given 
relatively little attention. Price /10/ perhaps best sums up the 
prevailing feeling when he says: 
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This is a field of work in which the theory was 
largely established well before any experimental 
verification was possible 5 and in which specta- 
cular discoveries were neither expected nor 
obtained. 

In 1962 5 however , Alekseyeva efc al fill reported such a 
discovery. They announced that their theory predicts instead of 
a flat plateau, a decrease of the ionization loss at very high 
velocities. Furthermore , they stated that their experimental work 
was in agreement with their theory. Their data show a drop of 
several per cent in the blob density of electrons in the region 
between ~ 200 and — 600 ~ *’) m In work done 

at the same time, the data of Stiller / 28 / show a slight tendency 
to "peak" at 'P — 750. Since previous data above — 200 

had a statistical accuracy no better than 2?4, earlier experiments 
w 7 ere able neither to confirm nor to refute the Russian theory. 

With these facts in mind, we decided to attempt the task of gather- 
ing sufficient data to provide a statistical accuracy of 1% at 
high values of nr . Since the area of high had pre- 

viously been neglected, we decided to concentrate our efforts in 
this area and to take only enough data below 'P - 100 to 
establish the fact that our measurements were sensitive enough to 
detect any variation in the curve. The precise details of the 
formulation of the experiment are given in the next section. 

We have multiple scattered a total electron track length of 
84.45 centimeters and counted 191^195 electron blobs. We find 
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no significant evidence for a departure from a flat "plateau" in 



the region 



T > 100. 
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2 . 



Formulation of the Experiment 

After the purpose of this experiment had been decided, it 



was necessary to formulate our exact procedures* Since we wanted 



pair production, using the gamma rays from the decay of neutral 
pions for initiation of the pair production* 



The neutral pions are obtained from interactions cf high 
energy primary beam particles with nuclei in the emulsion* 

K. I. Alekseyeva et a'l / 2 7/ used this same procedure, using 19.6 
Rev. protons as the primary beam particles* 

The choice of primary beam particles was made with two 
criteria in mind. First, the beam energy should be as high as 
possible, since the number of neutral pions formed in the primary 
star (and hence the number of pair produced electrons) is a func- 
tion of the beam energy* Second, Shapiro / 29 / indicates that the 
particles used for normalization should have HT > 100. To satisfy 
these two conditions, pellicles were used that had an incoming 
HT beam of 16.2 + 0.6 Eev. This corresponds to ^ 115 
and enabled us to use the primary beam particles as our normaliza- 
tion particles. Because of the source of electrons, an area one 
centimeter square located seven centimeters from the primary 



to measure particles with • as high as possible, we chose to 
work with electrons. One source of these high T electrons is 



-77 




T 
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beam entrance end was chosen. This area was then scanned for 
electron pairs. There were 42 electron pairs found in pellicle 
fj “132, 67 in pellicle *133, and 72 ir„ pellicle H *134. 

The pellicles we used were from a 600 micron Ilford K-5 
emulsion stack which was exposed to a 16.2 + fc 6Bev beam 

at GEM in 1960. This stack was flown to Berkeley, California, 
and developed there four days after exposure. Of the 180 pelli- 
cles in the stack, numbers 121 through 180 are at the United 
States Naval Postgraduate School on loan from Dr. Walter H. Bark as 
at Lawrence Radiation Laboratory. 

Although it would have been desirable to measure the rate of 
ionization loss over the entire energy range with a single species 
of particle, electrons of energy less than 30 Mev are difficult 
to measure. This difficulty arises largely because at these energies, 
electrons suffer large multiple scatterings. Because of the strong 
scattering: (1) it is difficult to contain the electron in an 

emulsion for any significant length; (2) the probability of in- 
advertantly shifting to another track is high; (3) it is difficult 
to separate the background count from the true count; and (4) the 
measurement of the electron track length becomes somewhat sub- 
jective. We therefore used pions and protons from the primary 
stars as our source of data for HT K 40. Although it was not 
our intention to measure accurately the details of the rise to 
plateau, we wanted to be sure that our measurements were sensitive 
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enough to detect a rise. 

In order to calculate the pv values of the electrons, protons, 
and pions, a program was written for the CDG 1604 computer at the 
USNPGS. A discussion of this program, along with the program flow 
chart and printout are included in the appendices 0 Three physi- 
cists did the multiple scattering on each track, scattering each 
electron track at least twice to prevent misidentification. For 
fast electrons the determination of pv is complicated by the radia- 
tive energy loss along the electron path. In order to minimise 
errors resulting from this effect, no single electron track was 
■scattered for more than 1.5 centimeters, and most were scattered 
for only 0.5 centimeters. Each computer program output was 
analyzed for bremsstrahlung effects (see appendices), and if a 
bremsstrahlung was suspected, the track was discarded. Furthermore, 
the same segment of electron track used in the measurement of vv 
was used in the measurement of blob density. Thus if some un- 
detected bremsstrahlung did occur along the electron track, the 
blob density and the value of would be averaged ever the same 
range of pv. 

As a measure of the rate of ionisation loss, we have used 
blob density rather than grain density. The greater information 
content of ether track parameters /22, 26/ is offset for near- 
minimum tracks by the ease and accuracy with which blob density 
data can be gathered. Accuracy in this sense means reproduci- 
bility. Price /10/ has stated that it is possible to obtain 
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results by blob counting that vary by no more than 0.5% between 
observers. 

Although the rise from minimum to " plateau" might have been 
different had we used grain density, the shape of the "plateau 11 
should not change. That is, if a drop from "plateau" exists when 
grain density is used as a measure of ionization, it will also 
exist when blob density is used. Therefore, all data in this paper 
is presented in terms of blob density. 

Two scanners were given the task of blob counting all tracks 
used. To prevent the need for calibration between the two scanners 
each had two pellicles of her own. All blob counting in each pelli 
cle was thus done by one observer. The electron and pioa tracks 
were normalized to tracks of the primary pion beam at the same 
depth in the emulsion. Periodic recounting of random tracks showed 
that several months were required to reach a satisfactory reproduei 
bility of 1%. Because of this effect, much of the data gathered 
in the early part of this experiment were eventually discarded. 

The data presented in this paper are from three pellicles. Addi- 
tional data from the fourth pellicle will be published at a later 
date . 
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3. Measurement of particle velocities 

Multiple scattering measurements were made on each track to 
estimate the product of momentum and velocity, pv. The velocity 
and the corresponding value for were then calculated. Scatter- 

ing measurements were performed on proton, electron, and pion 
tracks using a Koristka R-4 microscope. The calculations were 



accomplished on a CTC 1604 computer; a program description, flow 
chart, and listing for which are included in the appendices. 

The same basic computation was used to estimate the value of 
pv for all particles. The expression /3<0/$ 



pv - 



K 

CO 




573 0 



(in Me.) 



was used for calculations. The scattering constant, 9 is chat 
appropriate for a cut-off without replacement at four times the 
true mean second difference, D. The cell length, s, is in microns. 

The scattering constant was evaluated from the expression; 

k* = £7S + <2272 Lo<j w (Ss'Jj 

where the constants are those given by Scott /31 / and adjusted 
by Barkas /30/. The equivalent cell length, s\ corrects for 
the dependence of K _ on particle velocity and charge. The same 
references supply the expression for the equivalent cell length; 

s « (0.23 + 0.77 v 2 /c 2 )(s') 
for singly charged particles. 
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Noise was removed from the scattering measurements for each 
track by two distinct procedures. The methods used were; (1) cell 
length variation /32,3Q/ using cell lengths of one, two, and three 
times the prime cell lengthy and (2) subtraction of a constant 
noise appropriate to each observer. 

Using cell length variation the nolss can be estimated from; 



Noise squared (£,j) 




J 

values of the measured second differences using ceil lenghts of 
i and j times s. The indices i and j vary from one to three; 
however x J • The factor of 0.14 in the exponent comes 

from considering the scattering constant as a function cf cell 
length. The noise was calculated using cell overlap combinations 
of (1,3), (1,2), and (2,3) times the prime cell length. The 
three values were compared for consistency but the value of noise 
using cell overlap of (1,2) was the only one used in estimating 
the true signal. 

Repeated scattering measurements were performed on the 16 
Rev primary picms by all observers. For picas of this energy, 
apparent scattering is almost completely noise. Using this 
noise, q, the true mean second difference for an electron track 



are the means of the absolute 



In this expression ^2 Ctnd 
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was estimated from: 




where ^ is the measured total mean second difference. Prime 
cell lengths of 100 s 250 5 and 500 microns were used with the 
resultant noise found to be nearly independent of these cell 
lengths . 

All tracks were scattered twice , usually by different observer 
Tracks which did not yield results in statistical agreement , 
either between observers or method of noise removal , were dis- 
carded or rescattered. 

For fast electrons s the determination of pv Is complicated 
by the radiative energy loss along the path of the electron. To 
minimize this, no elec cron track was scattered more than 1.5 cm. 
Further, any electron track which showed evidence of a change in 
pv along the scattered length, either as a noticeable single 
scatter or as indicated by the computer program, was discarded. 

The computer divides each track into segments ten cell lengths 
long to facilitate the detection of changes in pv possibly caused 
by bremsstrahlung. (See app, I) 

Of a total of 362 electron tracks (181 pairs) located in 
the scanned areas of the three plates, 225 were discarded because 
of detectable bremsstrahlung or because the length of track in 
one pellicle did not allow a minimum of 20 prime cell lengths. 

In addition to scattering each track twice, requiring agree- 
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ment between methods in determination of B, and requiring that no 
detectable bremsstrahlung occur in the used portion of the track; 
all tracks used met the conditions shown in Fig 0 1. 





Min. 

sig/n 


Min No, 
of S 


S used 
(micron) 


Optimum 
No, of S 


pv' ^ 2, 00 


1.5 


20 


iOO 


40 


200 < pv < 1000 


1.5 


20 


250 


40 


pv > (000 




20 


500 


30 



Fig. 1 Useability Criteria 



Fig. 3 / 33 / was used to estimate the minimum cell length to use for 
scattering to expect a signal to noise ratio ^ 1.5 for a 

specific value of pv. 



The standard deviation, Q* , of pv determinations was esti- 



mated using CP s fiFpn /32/; where n is the number of prime cells. 

r “ 

/\ is the degree of overlap, and f is taken from Fig. 2. 



A 


J 


2 


3 


f 


1 


2.022 


2.026 



Fig. 2 Overlap Form Factors 



Using the criteria in this section, the estimated errors in pv 



ranged from 15 to 30 percent. 
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FIG. 3 



4. 



Measurement of ionization 



Blob density has been used as the measure of ionization in 
this experiment. The greater information content of other track 
parameters, such as grain density or mean gap length /34/ s Is 
offset for near minimum tracks by the ea.se and accuracy with which 
blob density data can be gathered. Bausch and Lomb miscroscopes , 
with magnification of about 16 x 100, equipped with coordinate 
stages by the Lawrence Radiation Laboratory at Berkeley, California, 
were used for blob counting. 

The energy of fast electrons decreases exponentially with 
increasing track length in the emulsion. As the blob density is 
also a function of energy, care was taken to assure the scattering 
measurements for second differences and the average blob densities 
were obtained over the same portion and lengch of each electron 
track. Therefore if some undetected radiative energy loss occurred 
along an electron track, the ionization loss and the value of T* 
are averaged over the same range of pv. 

Blob count calibration curves for the three plates used are 
shown in Fig. 4. These curves describe the variation of blob 
density with depth in the emulsion for each of the pellicles. To 
obtain these plots, the pellicle thickness was divided into tenths 
and then the blob density of che primary pious within each tenth 
of pellicle thickness was determined. Each point represents 
about 12 tracks or 5,000 blobs. These tracks were all located 



in the same region as the electron pairs. The variation in these 
curves for pellicles processed together s clearly shews the need 
for careful calibration of each pellicle. (See subhead 7 5 Critique) 

All blob density data were collected by twe scanners. All 
counting in a given plate was done by one. Bath scanners were, 
required to return periodically to their "normalization tracks'” 
(primary pions) to check for any change in subjective criteria used 
in blob counting. Changes did occur s and much of the data gathered 
in the early part of the experiment had to be discarded because 
of this. Several months were required to reach reproducibility of 
about one percent. It should be noted the two scanners had a 
continued systematic difference of about 3 % in their counts of the 
same tracks. As a consistency check the shape of one of the curves 
of Fig. 4, that for £}-134, was reproduced by each scanner. 

The electron and pion blob densities were normalised by form- 
ing the ratio of blob density of a track to that of the primary 
pions at the same depth in the emulsion. To select the primary 
pion blob density to be used for normalizing a given track 9 the 
initial and final depths in the emulsion for that portion of the 
track used in data gathering were determined. The calibration 
curves were considered as linear sections and a weighted average 
track depth was used to establish the depth in the emulsion of the 
primary pion to be used for normalization. 
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1 7 







5. 



Results 



The effect which this experiment was designed to detect was 
a 5% variation in ionization loss at values of T greater than 
200. For T" greater than 10, data were processed to provide a 
statistical uncertainty on all points plotted in Fig. 12, 13, and 
14 of 2% or less. At lower r 'T values, no statistical uncertainty 
exceeds 3%. 

To provide this degree of certainty, several tracks were 
combined to produce each data point, but each point represents 
only one type of particle. Tracks were added together, beginning 
with the lowest value of T , until the sum of the blob count, 

N, for each point was such that 100 x/ % NJ /£? N = 2% or 3% as des- 
cribed above. 

Since the 16.2 BeV pions ( Hf ^ 115) to which the stack was 
exposed have an ionization loss which differs by less than 1% 
from the ionization loss of the plateau, the blob count for these 
tracks was selected as the value to which all others should be 
normalized. The effect of normalization on the data for each plate 
is illustrated by comparing Figs. 6, 7, and 8 with Figs. 9, 10, 
and 11 for plate numbers 132, 133, and 134 respectively. 

The following relations describe how the normalized mean 
blob count, c, the ordinate of Figs. 12, 13, and 14, was 
determined for each point. 
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N = Number of blobs in the ith track 
1 

N = Number of blobs in the same length of 16.2 BeV pion track 
o 

nu ■ Number of blobs per 100 microns of track 

n = Number of blobs per 100 microns of 16.2 BeV pion track 
o 



n. = N./m. 
i li 



c. = n./n = N./m.n or c, = J ~ ^ assuming N. has 

i i o 110 i nrr n i 



N i t 



1 o 



a statistical uncertainty of 



(J- 2 (c. ) ■ N./m 2 n 2 

U N 1 1 1 o 

The c is the average of the c^'s weighted according to the inverse 

square of the statistical uncertainty. The summation is over the 
number of tracks involved . j 2 1 






r 



T z (Ci) 



urn 




— 



TciNj. 



2 2 

since ( n 0 / n £ ) adjusted only the third 



E Ni 

decimal place of the normalized blob count. Therefore each track 
is weighted according to the number of blobs counted. 

The value of HT for a given point represents the arith- 
metic average of the values found for each track represented by 
that point. On the plateau, the maximum spread of the values 
represented by any given point is 22% of the value at that point. 



19 



(Fig. 14, plate 134, = 2600.) A typical spread is about 107 o . 

Thus the points are separated by a distance somewhat greater than 
the error in T . 

These processes were used for determining the T and normal- 
ized blob count of each point for each of the three plates, sepa- 
rately, to produce the curves of Figs. 12, 13, and 14. Combining 
these plates by superimposing the plateau regions produced the 
curve of Fig. 5. This method of combination was dictated by the 
overwhelming statistical weight associated with the information 
in the plateau region as compared to the information in the re- 
mainder of the curve. (For electrons, a total track length of about 
84 cm. was examined.) 
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Fig 
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Fig. 10 






6. Comparison with theory and other experiments 

As was described in Subhead 1, the relativistic rise to 
plateau has been investigated by many persons. In 1950 9 Pickup 
and Voyvodic /l 1 / reported indications of a relativistic increase 
in grain density of about 10% to a plateau starting at HP — 20, 

Similar results were obtained by Daniel et al /15/. In 1953 
Stiller and Shapiro / 1 7/ reported a 14 + 3 % rise to plateau. In 
our experiment, as in Stiller and Shapiro's, blob densities were 
substituted for grain densities to simplify counting and facilitate 
reproducibility among observers. The magnitude of the relativistic 
rise of blob count determined was 14 4- 2%. In our experiment . 

The rate of rise to the plateau and the constant value on 
plateau indicated by our results is essentially the same as that 
determined by Stiller and Shapiro and many others but extends the 
investigated region on the plateau to higher values of Hf . The 
effect noted by Alekseyeva et al /27/ was not observed. 

Our data are entirely consistent with the theory as determined 
by Halpern-Hall / 7/ , and Sternheimer /9/ s with, the plateau begin- 
ning at T greater than 100. This indicates a deviation from 
the theories of Daniel et al / 15 / and Morrish /14 / which predict 
a plateau beginning about T = 20. 
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7. Critique 

Although the data seem to be in good agreement with currently 
accepted theory, there are possible sources of error which could 
obscure a departure from a flat plateau of ionization loss* These 
are the determination of from multiple scattering measure- 

ments and the variation of blob density with depth in the emulsion 0 

To check for systematic errors in our scattering data*, the 
pellicles were taken to the Lawrence Radiation Laboratory at 
Berkeley and a random sample of tracks was scattered on the Koristka 
MS-2 belonging to the Barkas group. Consistent agreement within 
statistical error was found, and we consider this the best evidence 
we have for thinking that our estimates of *V are not systemati- 
cally in error. 

Since our data were insufficient to allow restricting all 
tracks to a small fraction of the emulsion thickness , a correction 
was required to compensate for the variation of blob density 
with depth. (See subhead 4.) Although pains were taken in the 
normalization, one must recognize that the variation with depth 
is a severe handicap. Some encouragement comes, however, from 
the consistent results obtained from pellicles with rather 
different normalizing curves. 
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APPENDIX I 



SCATTERING PROGRAM 

This program was written for the purpose of calculating the 
pv of a singly charged particle in a nuclear emulsion. The method 
of multiple scattering is described in the section on determina- 
tion of velocity. Our program was written in Fortran 60 for 
specific use on a Control Data Corporation 1604 computer. A dis- 
advantage for use as a general scattering program is that in order 
to compute an equivalent cell length, an estimate of v/c must be 
made. The data printout format was designed so that the actual 
computer paper could be filed and kept as a record. During the 
course of the experiment, printouts were added when deemed advan- 
tageous; until in its final form, the program outputs sufficient 
data to analyze the run. 

The following are the program inputs. On the first card of 
each run: 

Columns 1-3 Pellicle number 
Columns 5-7 Event number 
Columns 9-10 Prong number 
^Column 11 Scanner identification 
Columns 12-17 Date of scattering 
Columns 18-22 Primary cell length 
Columns 23-25 Number of Y coordinate readings 
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Columns 26-30 Estimate of v/c 



Columns 31-35 Microscope calibration 

* This number is used for filing and for choosing a personal 
noise reading written into the program. 

On the second card of each run are the following: 

Column 1 Number one if final track for this computer 
run; otherwise number zero. 

Column 2 Number one if computer is to calculate pv 
with both calculated and personal noise; 
otherwise number zero. 

Column 3 Always zero. 

On the following data cards there are eight coordinate readings to 
each card, one reading every ten columns. 

After reading the input cards, the computer calculates a noise 
reading; and using this noise reading, a value of pv. After a 
first run the program selects the personal noise reading for the 
individual doing the scattering, returns to the beginning, and re- 
calculates pv using this personal noise reading. The personal 
scanner noise was determined by having each person scatter the 
16.2 Bev pions of the primary beam. Cell lengths of 100, 250, and 
500 microns were used to scatter the pions; and since at these 
cell lengths the signal from a 16.2 Bev particle is virtually non- 
existent, any reading is noise. Enough tracks were scattered by 
each person to determine a mean value of individual noise. The 
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noise determined by this method was found to be independent of the 
cell length used. 

The program outputs were selected so that crosschecks could 
be made. The outputs were: 

1. The number of Y coordinate readings. 

2. The Y coordinate readings. 

3. All calculated second and third dif ferences--if a 
second or third difference was cast out by the guillo- 
tine method, the corresponding second or third dif- 
ference will print out as zero. 

4. Number of readings cast out by the guillotine method -- 
This output, together with output 3, can sometimes 
pinpoint a bremsstrahlung along an electron track. 

5. Mean second and third dif ferences--The calculation of 
second and third differences is made for primary cell 
lengths, double overlapping cell lengths, and triple 
overlapping cell lengths. There are six corresponding 
values of mean differences and cast outs. 

6. The three calculated values of noise—Noise was calcu- 
lated using primary and double cell lengths, primary 
and triple cell lengths, and double and triple cell 
lengths. Only the first value was used to compute pv, 
but the other two were printed out as a consistency 
check. 
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7. The actual second difference signal--This value is the 
square root of the mean second difference squared using 
primary cell lengths minus the calculated noise squared 
using primary and double cell lengths. 

8. The signal using the mean third difference and the 
ratio of outputs 7 and 8. 

9. The three values of pv (using the three mean second 
differences), the estimated statistical error, and the 
signal to noise ratio, all using calculated noise. 

10. The three values of pv, the estimated statistical 
error, and the signal to noise ratio, all using personal 
scanner noise--If outputs 9 and 10 did not agree within 
the estimated statistical error, the track was either 
rescattered or discarded. 

11, The mean second differences for primary cell lengths 
divided into segments of ten consecutive second dif- 
ferences --This output was found to give the best 
indication of a bremsstrahlung along an electron 
track. A bremsstrahlung can be suspected when the 
mean second difference increases to a higher value 
between successive segments, and then remains at 
this higher value in succeeding segments. 

There is one noteworthy aspect to this program. In some 
situations the calculated noise can be a negative number. Since 
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V 

this situation has no physical significance, the calculated noise 
is set equal to zero. When the computer is later asked for a 
ratio of signal to noise, it can be troublesome. This situation 
can be programmed out, or it can be handled as we did. When the 
calculated noise from primary and double cell lengths prints out 
as negative, the signal to calculated noise ratio is ignored-- 
it is always infinite. 



37 



APPENDIX II 



IOGRAM 







LOW Lhart 




38 




39 









V 




40 




41 







42 





-o 

43 





44 





4S 



APPENDIX III 
PROGRAM LISTING 



PROGRAM SCAT 

PROGRAM SCAT SYMBOL DEFINITIONS 

IDEV FIRS! 3 i PELL NBR 2ND 3, STAR NBR, LAST 2, PRONG NBR 
IOSC SCANNER NUMBER 

1 0 ATE MONTH, DAY, YEAR DATA TAKEN 
S CELL LENGTH IN MICRCNS 

N NUMBER OF READINGS IN EV.NT 

BETA ESTIMATED RCTA Al INPUT 

LAST NON-ZERO FOR LAST EVENT IN TUN 
Y ( I ) ORDINATES IN MICROMETER UNITS 
D I F F 2 ( M , I ) IS I-TH 2ND DIFF F' R CELL LENGTH MS 
DIFF3IM.I) IS l-TH 3RD DIFF FOR CELL LENGTH MS 
DBAR2IM) IS MEAN 2ND DIFF FOR CELLS MS, WITH CUT-OTF 

DBARi(M) IS MEAN 3RD CIFF FOR CELLS MS, WITHOUT CUT-OFF 

SIGSQ(N) IS SQUARE GF NOISE CORRECTED SEC DIFF FOR CELL NS, 

USING NOISE FROM SINGLE AND DOUBLE CELLS 
ESS1 IS EFFECTIVE S FOR CALC OF KCO 

FCOINI IS KCO FOR CELL NS, IN FLT PT 
KHEK(N) +1 IF PRC CALC FOR CELL NS 

0 IF NG BECAUSE NOISE GREATER THAN SIGNAL 
-1 IF NG BECAUSE SIG LESS THAN 2 X NOISE 
KHEK3 +1 IF 3RD DIFF CALC OK, 0 IF REJECTED 

46 PRINT 47 

47 FORMAT (9H1 J.N. DYER //) 

DIMENSION Y ( 200 ) ,KNIX< 6 ) ,DIFF2< 3, 19 3 ) , D I FF3 ( 3 , 1 98 ) , DBAR 2 ( 3 ) , 

1 DB AR 3 ( i ) ,ERRSQ(3,3»,SIGSG(3),FCC(3) , KHEKI 3 ) , RDBAR ( 6 ) , RE RRSQ ( 3 , 3 ) , 
2SIGMIC(3),PBC(3),RBARS2(3),RATIO{ 31 , 

3SEG ( 20 ) 

READ IN DATA 

1 FORMAT (110, II, 16, F5.C, 13, FG.3, F5.3) 

2 FORMAT ( 8F1 C . 0 ) 

3 FORMAT (11,11,11) 

READ 1, IDEV, I DSC , IDATE, S, N, BETA, CALIB 
READ 3, LAST, INKER. ITEMP 
RE AO 2, ( Y I I 1 , I = 1 , N ) 

REARRANGE FOR POSSIBLE DUPLICATES 
11=0 
N1 =N 

DELT A=0 . 

00 3000 I = 1 » N 

IF I Y ( I ) ) 3002,3001,3001 

3001 Y ( I ) =Y( I >+ OELTA 
Y ( I - I 1 I = Y( I) 

GO TO 3000 

3002 N 1 = N 1 - 1 
11= 11+1 

DELTA = Y( I) + Y I I — 1 ) 

3000 CONTINUE 
N=N 1 

CALCULATE SECOND DIFFERENCES 
00 4 M = 1 , 3 
NSTOP = N-2*M 
DO 5 1=1, NSTOP 

L = I + 2 * M 
K = I +M 

DIFF2(M,I) = Y(L)-2.*YCK)+Y( I ) 

5 CONTINUE 

4 CONTINUE 
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CALCULATE MEAN SECOND OIF FERtNCES WITH 
CUTOFF AT 4 DbAR 
PO 6 M=1,3 

TOTAL = C. 

NIX = 0 
KNIX(M) = 0 
NS TOP = N-2»M 
102 DO 7 1= UNSTOP 

TOTAL = TOTAL + A B SF ( D I FF2 ( M , I > ) 

7 CONTINUE 

A = NS T OP - KNIX(M) 

DRAR2IMJ » TOTAL/A 
RDBAR(M)=D3AR2( M)*CAL IB 
DO Q 1= UNSTOP 

IF (4.*0BAR2(M)-A8SF(0I FF2(M, I>) ) 9,9,8 
9 DIFF2IM, I ) = 0. 

NIX = NIX+1 
KNIX(M) = KNIXIM1+1 

8 CONTINUE 

100 IF (NIX) 6,6,101 

101 NIX = 0 
TOTAL = 0. 

GO TO 102 

6 CONTINUE 

CALC 3RD OIFFS 
DO 11 M=1.3 
NS TOP = N-3»M 
DO 10 1= UNSTOP 

01 FF3 ( M» l ) = Y ( I + 3*M ) - 5 . » Y ( I +2* M ) +3 .* Y ( I +M )- Y ( I ) 

10 CONTINUE 

11 CONTINUE 

CALCULATE MEAN THIRD DIFFERENCES WITH 
CUTOFF AT 4 DBAR 

DO 12 M=1,3 
TOT AL = 0 . 

N I X=0 

KNIX ( M ♦ 3 ) = 0 
NS TOP= N- 3 * M 
17 DO 13 1=1 , NSTOP 

TOT AL= TOTAL + ABS F ( D I F F 3 ( M , I > ) 

13 CONI INUL 

B= NSTOP- KNIXIM+3) 

D B A R 3 ( M ) = TOTAL/R 
RDDAR(M+3)=DBAR3(M)*CAL IB 
DO 14 1=1, NSTOP 

IF (4,*DPAR3(M)-ABSF( DI FF3(M, 1 1 ) ) 15,15,14 
'15 DIf F?(M, I 
N I X = N I X + 1 

KNlX(M+3)=KNIX(M+3)+l 

14 CONTINUE 

IF (NIX) 12,12,16 
16 NIX =0 
TOTAL =0. 

GO TO 17 

12 CONTINUE 

CALCULATE NOISE SQUAREO 
A = 3 . + U /LOGF ( 5. *S ) 

DO 20 M = 1 ,2 
DO 21 L = 2,3 
IF (M-L) 24,21,24 
24 FN = L 
FM = M 

ERRSQ(M,L)=DBAR2(M)«*2.-( ( CBAR2 ( L )* * 2 .-DBAR2 ( M ) **2. 1 / ( ( FN/FM 1 *« A- 

1 1 . 1 ) 

RERRSQ(M f L) = (E«RSQ(M,L) )* ( CAL IB**2. ) 

21 CONTINUE 
20 CONTINUE 
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CALCULATE PBC FRCM 1 AND 2 CELLS 
USING NOISE FROM 1 AND 2 CELLS 
22 DO 55 1=1,5 

I F ( I TEMP I 122, 122,1122 
1122 ERRSQ ( 1 ,21 =PERSER 
GO TO 31 

E R R S 0 ( M , N ) IS NOISE SQUARED FROM CELLS MS ANO NS 
122 IF ( E R R S Q ( 1 » 2 ) ) 30,31,51 

30 FRRSC (1,2) = 0. 

31 IF ( DRAR2 (11**2. - ERRSC(1,2)) 32,32,35 

33 SIGSG(l) = DB AR2 ( I ) » * 2 • - ERRSQI1,2) 

GO TO ( 34 1 , 3*40 » 34 ) , I 

341 ESS = S 
GO TO 342 

340 ESS = 2 . *S 
GO TO 342 

34 ESS=3.*S 

342 ESS1 = ESS/ ( .23 + . 77*BETA»*2 . I 

FCC(I) =SQRTF( 675.*! . 09 + . 2 72*LOG 1 OF ( 5.*ESS1 ) ) ) 

SIGMIC(I) = CALie*SQRTF(SIGSQ(I) I 

PBCIII = ( FCO! I )«ESS**1 .51/(573. *SIGMIC( I) ) 

RAT I 01 1 ) = SQRTFISIGSQl Il/ERRSQI 1,2 ) ) 

KH EK ( II = +1 
GO TO 55 

32 KHEK ( II = 0 
55 CONTINUE 

IF(ITEMP) 2000,2000,2020 

CALCULATE SIGNAL FROM 3RD DIFFS 
>000 IF!UBAR3(2)**2-DRAR3( 11**2) 200,200,201 

200 KHEK3 = 0 
GO TO 202 

201 SIG3 = SCRTF! ( DBAR3! 2 ) * »2.-DRAR 3( 1)*«2. 1/11.5*1 (2. **A)-1. ))) 

SIG3=SIG3*CALIP 
KHLK3= +1 

202 GO TO (60,61,63,64,66,67,69,70), IDSC 

60 I SC AN=8H JOHN 
GO TO 72 

61 ISCAN=8HFRLD 
GO TO 72 

63 I SC AN = 8HGA I L 
GO TO 72 

64 ISCAN=8HSHEILA 
GO TO 72 

66 ISCAN=CHHARRY 
GO TO 72 

67 I SCAN=8HRUSS 
PEKSLR =438.0745 
GO TO 72 

69 I SC AN=8HSWEDE 
PERSER =654.403 
GO TO 72 

70 I SC AN = 8HD I CK 
PERSER = 654.408 

72 PRINT 73 

7 30F0 *MAT ( 1 OX , 1 OH EVENT NO. , 10X , 8HSCANNER , 10X,6H DATE, 10X,5H S 
1 1 0 X , 3H N) 

PRINT 74. IDEV, ISCAN, IDATE, S,N 
74 FORMAT ( lOX, ll0,10X,A8, 10X, I 6, 1 OX ,F 5 . 0, 10X, 13 //) 

PRINT 75 

750F0RM AT ( 15X, 5H Y ( I ) , 8 X , 8H 1 D I F2 ( I ) , 7X , 8H2D I F2 l I ) , 7X , 8H3DI F2 ( I ) , 8X 
18H10IF3( I )»7X,8H2DIF3( I ) » 7X, 8H3DI F3 ( I)) 

PLACE MATRIX IN HERE 
PRINT 76, Y( 1 ) 

76 FORMAT ( 10X, F10.0) 

PRINT 77 , Y ( 2 ) » D I FF2 ( 1.1) 

77 FORMAT! 1 OX, F 10.0, 5X. FlO. 0) 

PRINT 78, Y(3),DIFF2( 1 , 2 ) , C I F F2 ( 2 , 1 ),DIFF3(1.1 ) 

78 FORMAT ( 10X.F10.0,5X,F10.0,5X,F10.0, 20X,F10.0) 

i PR I NT 79 ,Y ( 4 ) , 01 FF2 ( 1 ,3),0IFF2(2,2),DIFF2(3,1),DIFF3(1,2),DIFF3I2, 

79 FORMAT! 1 OXF 1 0. 0 , 5XF 1 0 .0 , 5XF 1 0. 0 , 5XF 1 0. 0, 5XF 1 0. 0 , SX, F 1 0. 0 ) 
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K=N-9 

PRINT ('5 » ( Y ( l +4 ) , u l F F2 ( 1 , I+3)»D1FF2!2» I +2 ) » 0 1 F F 2 ( 3 , I + 1 1 , 

1 D l F F3 ( 1 , I+2),DIFF3(2* 1 + 1) , 0 IFF 3 (3,1 ) , 1 = 1,K ) 

85 FORMAT ( 10XF10.0,5XF10.0,5XF10.0,SX=10.0,5XF10.0,5XF10.0 ,5XF10.0) 
K= N-7 
L = N-6 

PRINT 62,!Y(J+3)»DIFF2( 1,J+2)»D1FF2(2,J+1)»DIFF2(3»J), 
ll)IFF5( 1 . J + l >,01FF3(2,J),J=K,L) 

62 FORMAT! 1 0X F 1 0 . 0 , 5X F 1 0 . 0 , 5XF 1 0 . 0 , 5XF l 0 . 0 , 5XF 1 0 . 0, 5XF 1 0.0 ) 

K = N-4 

PRINT 65, Y(K+2) ,DIFF2< 1 , K + 1 ) , D I F F2 ( 2 , K ) . D I FF3 ! 1 , K ) 

65 FORMAT! 1 0X F 1 0 . 0 , 5XF 1 0 . 0 , 5X F 1 0 . 0 , 20X F 1 0 . 0 ) 

L = N - 3 

PRINT 68, Y(L+2),DIFF2( 1 , L+ 1 ) , D i FF3 < 1 , L ) 

68 FORMAT! 1 OXF 1 0.0 , 5XF 1 0 .0 , 35XF1 0 . 0 ) 

PRINT 7 1 , Y ( N I 
71 FORMAT ( 1 OX , F 1 0 • 0 ) 

I A W AY = 7HC AS T CUT 

PRINT 80, I AWAY.KNIX! 1 ) , KNIX! 2 ) ,KNIX ( 3) ,KNIX (4) ,KN1X( 5 ),KNIX( 6) 

80 FORMAT! 1 OX , A7 , 1 4X , 1 5 , 1 2X , I 3 , 1 2X , I 3, 1 2X, I 3, 1 2X, 13, 12X, 13 ///) 

I S 1 G=6HS I GN AL 
IPBC=3HP8C 
I DB AR=4H0BAR 

OPR I NT 86 , I OB AR , 0RAR2 ( 1 ) ,DBAR2(2 ) , DB AR2 ( 3 ) , 0BAR3 ( 1 ) ,DBAR3(2 ) , 
10BAR3! 3) 

860F0RMAT ( 10X,A4,11X,F10.2,5X,F10.2,5X,F10.2,5X,F10.2,5X,F10.2,5X 
1F10.2 /) 

PRINT 35, RDBAR( 1 ) ,RD(?AR( 2 ) ,RDBAR ( 3) ,R08AR(4) ,RDBAR(5 ) , 

1 RDB AR ( 6 ) 

35 FORMAT ( 1 OX , 7HDBARM l C , 8 X, F 1 0 . 3, 5X ,F 1 0 . 3 , 5X , F 1 0 . 3, 5X , F 10 . 3 , SX 
1F10.3,5X,F10.3 /) 

I NO I SE = 8HN0 I SE SO 
I F IRST= 5HM 1 N2 
IStC= 5HM1 N3 
I TM I RD=5HM2 N3 
DO 43 M= 1 , 3 

RBARS2(M) = ( ( DBAR2(M)*CALIB)«*2. ) 

43 CONTINUE 
J0Y=8HDBARSQMI 

PRINT 44 , JOY ,RBARS2( 1 ) , RBARS2I 2 ) , RB ARS2 ( 3) 

44 FORMAT! 10X, A8,5X,F10. 5, 5X, F10.5,5X, F 10.5 ///) 

890PR I NT 87,INOISE, IFIRST, ERRSQ ( 1,21 , 1 SEC,ERRSQ! 1 ,3), I THIRD, 
1ERRSQ(2,3) 

87 FORMAT! 10X, A8, 5X, AS, 2X, F10.3, 5X,A5, 2X , F 1 0. 3 , 5X, A 5 ,2 X , FI 0.3 /) 
ICHECK=6HN0ISE 
JOY2=8HLRRSOMIC 

v OPRINT 45, J0Y2, IFIRST,RERRSO( 1 ,2 ) , IS EC , RERRSQ ( 1 , 3 ) , I TH1R D , 

1 RE RRSO ( 2 , 3 ) 

. 45 FORMAT! 1 OX , A8, 5 X , AS , 2X , F 1 0 . 4, 5X , A 5, 2 X , F 1 0. 4 , 5X , A 5, 2 X, FI 0. 4 ///) 

DO 4 1 M = 1 , 3 

IF ! KHL K ! M ) ) 97,98,97 

97 PRINT 3003,PBC!M).RATIO(M) 

3003 FORMAT! 1 OX , 3HPBC , 1 OX , FlO. 1,24H SIGNAL TO NOISE RATI0,F10.2 //) 
GO TO 4 1 

98 PRINT 42 

42 FORMAT ( 10X, 3HPBC, 10X,2PHN0ISE IS GREATER THAN SIGNAL //) 

41 CONTINUE 

LD I F = 3HBY 3DIFF 
IF ( KHEK3-0 )92»92»90 

90 PRINT 91 , I S I G ,LDIF,SIG3 

91 FORMAT! 1 OX , A6 , A8 , 6X, FlO. 3 //) 

GO TO 94 - 

92 I B A0=4HBARF 

PRINT 93, I S I G , LDI F, IBAD 

93 FORMAT! 1 OX , A6,.A8, 6X, A4 //) 

94 CONTINUE 

PRINT 120, SIGMIC(l) 

120 FORMAT ( 1 OX , 1 5HS I GNAL BY 2DIFF, 5X,F 1 0.3 // ) 

BUDGE = SIG3 / SIGMIC! 1 ) 

PRINT 121, BUDGE 

121 FORMAT ( 1 OX , 27HS I GNAL RATIO 3DIFF TO 2DIFF,F10.3 //) 
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CALCULATE OOAR FROM SEGMENTS OF 
10 D I F F2 ( 1 ) 

TOT AL=0. 

KA =0 
MA = N-2 
M I =0 

DO 110 I=>1,MA 

IF (OIFF 2 ( 1 , I) ) 111,112,111 

111 TOTAL = TOTAL ♦ ABSF ( D I FF2 ( 1 , II ) 

MI = MI +1 

IF ( MI - 10 ) 113,112,113 
113 IF (I - MAI 110,112,110 

112 IF (MI) 114, 110,114 
1 1 4 KA * KA + 1 

B = MI 

SEG (KA) = TOTAL/B 
MI = 0 
TOTAL = 0. 

110 CONTINUE 

CALC SIGNAL IN MICRONS FOR SEGMENTS 
DO 115 1=1 ,KA 

I F { ( SEG ( I ) )**2.-ERRSQ{ 1.2) ) 116,116.117 

117 SEG( I ) = (SQRTF( ( SCG( I ) )»«2.-ERRSQ( 1, 2) ) )* CALIB 
GO TO 119 

116 SEG (I 1=0. 

119 PRINT IIP, SEG ( I ) 

118 FORMAT 110X , 22HSEGMENT SIGNAL M ICRDNS , SX, F7 . 3 /) 

115 CONTINUE 

ESTIMATING ERROR IN P.'C 
2020 DIMENSION PBCERI3) 

PRI M=N 

DO 400 1=1,3 
TEMP = KM X ( I I 
FL = I 

F ACT 1 = 1. + .022MFL-1 . ) 

FACT2 = PRIM - (2.*FL+TEMP) 

DENOM = SQRTF( FACT 1 »F ACT2/FL ) 

PBCER(I) = PBC( I ) *1.37/CEN0M 
400 CONTINUE 

I F ( ITEMPI402, 402, 2030 

402 PRINT 403 

403 FORMAT! 10X»45HERROR IN MEV CALCULATED FROM PG264 CERN NOTES/) 
DO 406 J =1,3 

PRINT 404, PBCER(J) 

404 FORMAT! 10X, 17HPBC PLUS GR M I NUS 5X ,F 10.3/) 

406 CONTINUE 

I T L MP= 1 

IF (INKER) 18,10,22 

2030 PRINT 2031 

2031 FORMAT ( 1 OX , 24HPPC USING PERSONAL NOISE/) 

DO 2033 1=1,3 

PRINT 2032, PBC(l). PBCER(I), R AT 13 < I ) 

2032 FORMAT (10X, F5.0,10X, 13HPLUS OR MINUS 5X, F5.0, 

110X, 21HSIGNAL TO NOISE RATIO 5X, F5.2/) 

2033 CONTINUE 

18 IF (LAST) 48,46,48 
48 CONTINUE 
408 END 
END 
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